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Abstract 

The logarithmically enhanced ln(l/a) corrections to the para- and orthopositro- 
nium decay widths are calculated in the framework of dimensionally regularized non- 
relativistic quantum electrodynamics. In the case of parapositronium, the correction 
is negative, approximately doubles the effect of the leading logarithmic a^ln^(l/a) 
one, and is comparable to the nonlogarithmic 0{a'^) one. As for orthopositronium, 
the correction is positive and almost cancels the ln^(l/a) one. The uncertainties 
in the theoretical predictions for the decay widths are reduced. 
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Positronium (Ps), which is an electromagnetic bound state of the electron e~ and the 
positron e"*", is the lightest known atom. Since its theoretical description is not plagued by 
strong-interaction uncertainties, thanks to the smallness of the electron mass rrie relative 
to typical hadronic mass scales, its properties can be calculated perturbatively in quantum 
electrodynamics (QED), as an expansion in Sommerfeld's fine-structure constant a, with 
very high precision. Ps is thus a unique laboratory for testing the QED theory of weakly 
bound systems. 

The decay widths of the ^5*0 parapositronium (p-Ps) and ^5*1 orthopositronium (o- 
Ps) ground states to two and three photons, respectively, have been the subject of a vast 
number of theoretical and experimental investigations. The present theoretical knowledge 
may be summarized as: 
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are the lowest-order results. The 0{a) coefficients in Eqs. (|I|) |]T| and 
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read 
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Ao = -10.286 606(10). 
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The logarithmically enhanced a^ln(l/a) terms in Eqs. (|I]) and have been obtained in 
Refs. 0, ^, respectively. Recently, the nonlogarithmic 0(q;^) coefficients in Eqs. (|1|) 
and d^) P have been found to be 



5p = 5.14(30), 
Bo = 44.52(26). 



(5) 



Note that the light-by-light scattering diagrams have been omitted in Ref . . In the p-Ps 
calculation the diagrams of this type increase the coefficient Bp by 1.28, and, therefore, 
their contribution to the coefficient Bo is assumed to be relatively small. The p-Ps (o-Ps) 
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decays into four (five) photons, which are not included in Eq. lead to an increase of 
the coefficient Bp (Bo) by 0.274(1) (0.19(1)) 0. In 0{a^), only the leading logarithmic 
a^ln^(l/a) terms are known [Q. Including all the terms known so far, we obtain for the 
p-Ps and o-Ps total decay widths 

r*^ = 7989.512(13) fxs-\ (6) 
F*^ = 7.039943(10) /is-\ (7) 

where the given errors stem only from the coefficients Bp and Bo respectively and we 
postpone the discussion of total uncertainty of theoretical estimates to the end of the pa- 
per. The purpose of this letter is complete our knowledge of the logarithmically enhanced 
terms of 0{a^) by providing the coefficients Cp and Co, in analytic form. We also give 
order-of-magnitude estimates of the unknown coefficients Dp and Do- 
On the experimental side, the present situation is not entirely clear. Recently, the 
Ann Arbor group measured the p-Ps width to be 

r^^P = 7990.9(1.7) fis-\ (8) 

which agrees with Eq. (|]) within the experimental error. However, in the case of o-Ps, 



their measurements [10, 11 



r^^P(gas) = 7.0514(14) fxs~\ 
r^^P (vacuum) = 7.0482(16) /is-\ (9) 

exceed Eq. (|^) by 8 and 5 experimental standard deviations, respectively. This apparent 
contradiction is known as the o-Ps lifetime puzzle. On the other hand, the Tokyo group 



found PI 



r^"P(Si02) = 7.0398(29) fis-\ (10) 

which agrees with Eq. (|^ within the experimental error. Leaving this aside, the o- 
Ps results from Ann Arbor could be considered as a signal of new physics beyond the 
standard model. However, a large number of exotic decay modes have already been ruled 



out |jT3|. No conclusion on the o-Ps lifetime puzzle can be drawn until the experimental 
precision increases and the data become unambiguous. 

On the theoretical side, it is an urgent matter to improve the predictions of the Ps 
lifetimes as much as possible. Thus, one is faced with the task of analyzing the 0{a^) 
corrections, which is extremely difficult, especially for o-Ps. However, there is a special 
subclass of the 0{a^) corrections which can be analyzed separately, namely those which 
are enhanced by powers of ln(l/a) ?a 5. They may reasonably be expected to provide 
an essential part of the 0{a^) corrections. This may be substantiated by considering 
Eqs. (|1|) and (0) in 0(a;^), where logarithmic terms enter for the first time. In the case 
of p-Ps, 98% of the O(a^) correction stems from the logarithmic term. In the case of 
o-Ps the logarithmic term is not so dominant but still gives about 1/4 of the total O(a^) 



3 



correction. The origin of the logarithmic corrections is the presence of several scales in the 
bound-state problem. The dynamics of the nonrelativistic (NR) e'^e~ pair near threshold 
involves four different scales \Ti\: (i) the hard scale (energy and momentum scale like 



me); (ii) the soft scale (energy and momentum scale like /3me); (iii) the potential scale 
(energy scales like /3^me, while momentum scales like /3me); and (iv) the ultrasoft (US) 
scale (energy and momentum scale like /3^me). Here (3 denotes the electron velocity in 
the center-of-mass frame. The logarithmic integration over a loop momentum between 
different scales yields a power of ln(l//3). Since Ps is approximately a Coulomb system, 
we have P oc a. This explains the appearance of powers of ln(l/a) in Eqs. ([l|) and (H). 
The leading logarithmic corrections may be obtained straightforwardly by identifying the 
regions of logarithmic integration j^, ^, § . The calculation of the subleading logarithms is 
much more involved because certain loop integrations must be performed exactly beyond 
the logarithmic accuracy. 

In the following, we briefly outline the main features of our analysis. We work in NR 
QED (NRQED) [|15|, which is the effective field theory that emerges by expanding the 
QED Lagrangian in (3 and integrating out the hard modes. If we also integrate out the soft 
modes and the potential photons, we arrive at the effective theory of potential NRQED 
(pNRQED) [l^, which contains potential electrons and US photons as active particles. 
Thus, the dynamics of the NR e~^e~ pair is governed by the effective Schrodinger equation 
and by its multipole interaction with the US photons. The corrections from harder scales 
are contained in the higher-dimensional operators of the NR Hamiltonian, corresponding 
to an expansion in /3, and in the Wilson coefficients, which are expanded in a. In the 
process of scale separation, spurious infrared (IR) and ultraviolet (UV) divergences arise, 
which endow the operators in the NR Hamiltonian with anomalous dimensions. We use 
dimensional regularization (DR), with d = 4 — 2e space-time dimensions, to handle these 
divergences [|T6|, [1^ . This has the advantage that contributions from different scales are 
matched automatically. The logarithmic corrections are closely related to the anomalous 
dimensions and can be found by analyzing the divergences of the NR effective theory. In 
this way, we have obtained the leading logarithmic third-order corrections to the energy 
levels and wave functions at the origin of heavy quark-antiquark bound states , which 
includes the QED result p, ^, |] as a special case. Here, we extend this approach to the 
subleading logarithms in QED. Note that the NRQED approach, endowed with an explicit 
momentum cutoff and a fictitious photon mass to regulate the UV and IR divergences, 
has also been applied to find the third-order correction, including subleading logarithms. 



to the hyperfine splitting in muonium |20 



The annihilation of Ps is the hard process which gives rise to imaginary parts in the 



local operators of the NR Hamiltonian |^ . The decay width can be obtained by averaging 
these operators over the bound-state wave function. The hard-scale corrections, which 
require fully relativistic QED calculations and are most difficult to find, do not depend 
on /3 and do not lead to logarithmic contributions by themselves. However, they can 
interfere with the logarithmic corrections from the softer scales. Thus, the only results 
from relativistic perturbation theory that enter our analysis are (i) the one-loop hard 
renormalizations of the imaginary parts of the leading four-fermion operators, i.e., the 
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Born decay amplitudes, which are given by the coefficients Ap and Ao, and (ii) the hard 
parts of the one-loop 0(a/5^) operators p2|]. The missing ingredients can all be obtained 
in the NR approximation. These include (i) the correction to the Ps ground-state wave 
function at the origin due to the 0(a/5^) terms in the NR Hamiltonian, (ii) the 0(a/3^) 
and 0(a^/3) corrections to the leading four-fermion operators, and (iii) the correction due 
to the emission and absorption of US photons by the Ps bound state. 

The value of the ground-state {n = 1) wave function at the origin ipi{0) may be 
extracted from the NR Green function G(x, y, E), which satisfies the equation 



{He + An-E) ^(x, y , E) = ^^^'^ (x - y) 



where Tic is the Coulomb Hamiltonian and ATi stands for the terms of higher orders in 
a and j3. The solution of Eq. ([TT|) can be found in time-independent perturbation theory 
as an expansion in a around the leading-order Coulomb Green function. We thus obtain 

2 

the correction Aipf in the relationship |'?/'i(0)p = ipii^) (1 + Aipf), where V^f (0) is the 
ground-state wave function at the origin in the Coulomb approximation. As mentioned 
above, this analysis may be enormously simplified by the use of DR. Proceeding along the 
lines of Ref. |^, we thus recover with ease the well-known a^ln(l/a) terms in Eqs. (|I|) 
audi) ii. 



In ■ 
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(12) 



where S is the eigenvalue of the total-spin operator S. As mentioned above, A'lpl interferes 
with the one-loop hard renormalization of the Born amplitudes, Ap and Aq, to produce 
a^ln(l/a) terms. The resulting contributions to the coefficients Cp and Co read 2Ap and 
—Ao/3, respectively. 

The generic logarithmic 0{a^) correction A"iljf to |?/'i(0)p is generated by the following 
one- loop operators, given in the momentum representation with 0(e) accuracy. 
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(13) 
(14) 

(15) 



where 1/e = 1/e — 7£; + ln(47r), with ■je being Euler's constant, yU is the 't Hooft mass scale 
of DR, q is the three- momentum transfer, and E^' = —a^me/^ is the Coulomb ground- 
state energy. Equations (p!3[) and ([T^) give the hard |22| and soft [|l^ 0{a(3'^) contributions 
to the NR Hamiltonian, respectively. The US contribution given in Eq. ([TB|) arises from the 
emission and absorption of an US photon, which converts the on-shell Ps ground state into 
some off-shell state of the Coulomb spectrum, with three- momentum p, before it decays. 
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It is the only US contribution which can be represented by an operator of instantaneous 
interaction and thus give rise to logarithmic corrections. It has been found with the help 
of the method developed for the more complicated case of quantum chromodynamics in 



Ref. p3[, where it was applied to the on-shell renormalization of the heavy-quarkonium 
wave function at the origin. The singularities of the operators in Eqs. (p!^)-(p!5|) yield the 
logarithmic corrections which we are interested in. Up to their logarithmic dependences 
on and p^, these operators are of the (5-function type in coordinate space and, therefore, 



lead to additional singularities in the Coulomb Green function at the origin |]T9|. As a 
consequence, in the evaluation of the Green function in time-independent perturbation 
theory from Eq. ([Tl|) with AiJi,, AgTi, and A^sTi., overlapping logarithmic divergences 
appear in the part of the first-order term which corresponds to the interference of the one- 
photon contribution to the Coulomb Green function and the first terms of Eqs. (]T3|)-(|T5|). 
This results in the double-logarithmic contribution, which can be directly extracted from 
the coefficient of the leading double-pole singularity Since we are interested in the 



single-logarithmic contribution, we also have to keep the subleading terms in this analysis. 
The logarithmic corrections which are generated by the non-overlapping singularities can 
be obtained by putting fi = rrif. in the Coulomb Green function at the origin and = 
= — me-Ef in Eqs. ([T3|)-([T5D, and proceeding as in the evaluation of Eq. (|T^) . We thus 
obtain 



vr L 2 a 



184 2, ^ 
+ -ln2 



45 3 



+ (^y + ln2)s(S + l)]lni|. (16) 

The first term herein agrees with the corresponding terms in Eqs. and (0) p|, while 
the second one represents a new result. 

Another source of Q;^ln(l/a) terms is the 0(a/3^) corrections to the leading four- 
fermion operators. Since they do not involve the singular Coulomb Green function at the 
origin, there are no overlapping divergences, and we may simply read off the resulting 
a;^ln(l/Q;) terms from the poles of their US parts, which are given by the operator 

1 2a p^ + p'\ 
--^ 2 — ^4P,P,S. 17 

Here V4(p, p', S) is the local four-fermion operator which generates the leading-order decay 
widths. Taking the expectation value of Eq. (|l^) w.r.t. the ground-state wave function, 
one encounters power-divergent integrals |^|. They can be consistently treated within 



DR 0. This leads to the substitution p , p' — > meE^ in the matrix element. The UV- 



pole contribution of Eq. ([T7| ) is then canceled by the IR pole of the hard contribution ||24 |. 
This implies that the logarithmic integration ranges from the US scale a^uif. up to the 
hard scale me, so that 1/e should be replaced by 41n(l/a) [|l^. The resulting a^ln(l/a) 
corrections to the decay widths are spin-independent and read 

Ar„„^ro;.i^.„i. (is) 
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The last source of a^ln(l/Q;) terms is the 0{a^(3) corrections to the leading four- 
fermion operators. These corrections are non-analytic in and of the form 



7 /„ 32 . IpI + |p'| 



rrip 



^4(p,p',S), (19) 



where the first and second terms contained within the square brackets are the soft and US 
contributions, respectively. Ahhough the coefficient of V4(p,p',S) in Eq. (|1^) is finite, 
the matrix element of Eq. ([T9|) between Ps bound-state wave functions is logarithmically 
divergent in the UV region. In DR, the divergent part is given by the matrix element 
of Eq. ( p^D with |p|, |p'| replaced by mea/{ne). In contrast to Eq. (|17D , the logarithmic 
integration now ranges from the soft scale up to the hard one, and 1/e should be replaced 
by 21n(l/a;). The resulting contributions to the decay widths read 

Summing up the various a^ln(l/a) terms derived above, we obtain 

367 

= 2Ap + — - 10 In 2 ^ -7.919, 
90 

A 229 

a = -— + — -81n2^ 5.517. (21) 
3 30 ^ ' 

In the case of p-Ps, the new a^ln(l/a) term has the same sign as the a^ln^(l/a) one and 
exceeds the latter in magnitude. The sum of these two terms compensates approximately 
1/3 of the positive contribution from the nonlogarithmic 0{a^) term. As for o-Ps, the 
new a^ln(l/a) term cancels approximately 3/4 of the a^ln^(l/a) contribution. Our final 
predictions for the p-Ps and o-Ps total decay widths, including the multi-photon channels, 
read 

r*^ = 7989.620(13) /is"\ (22) 
r*^ = 7.039968(10) /is'S (23) 

which has to be compared with Eqs. (^ and (0). As before, Eq. (|2^) agrees with the Ann 
Arbor ^ measurement (H), and Eq. (|23| ) favours the Tokyo measurement (|13), while 



it significantly undershoots the Ann Arbor [[10| , |ri|] measurements 



The missing nonlogarithmic 0{a^) corrections in Eqs. (0) and (|^) receive contributions 
from three-loop QED diagrams with a considerable number of external lines, which are far 
beyond the reach of presently available computational techniques. In this sense, we expect 
Eqs. (^) and (p3|) to remain the best predictions for the forseeable future. However, we 
may speculate about the magnitudes of the coefficients Dp and Do- Two powers of a in 
these terms can be of NR origin. Each of them should be accompanied by the characteristic 
factor vr, which happens for the logarithmic terms. Thus, we estimate the coefficients Dp 
and Do to be a few units times vr^. This rule of thumb is in reasonable agreement with 
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the situation at O(a^), where we have Bp ^ 7r^/2 and Bo ~ 47r^. If the coefficients Dp 
and Do do not have magnitudes in excess of 100, then the uncertainties due the lack 
of their knowledge falls within the errors quoted in Eqs. ( p^ ) and (^). Then, our new 
results reduce the uncertainties in the predicted p-Ps decay width to 10~^yus~^. The main 
remaining theoretical uncertainty in o-Ps decay width is related to the unknown 0{a^) 
contribution of the light-by- light scattering diagrams which can be estimated as a few 
units times 10~^/is~^ on the basis of p-Ps result. Calculation of this contribution along 
with our present result will reduce the uncertainties in the predicted o-Ps decay width 
to 10~^/is~^. Further progress in our understanding the Ps lifetime problem crucially 
depends also on the reduction of the experimental errors, which now greatly exceed the 
theoretical ones. 

Finally, we note that the technique developed in this letter can also be applied to the 
calculation of the subleading logarithmic a^ln(l/a) terms for the Ps hyperfine splitting. 
This problem is of special interest because of the apparent discrepancy between the latest 
experimental data and the best theoretical predictions, which include the 0{a^) 



a 



corrections (see Ref. and references cited therein) and the leading logarithmic In 
term g. 
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Contract No. 05 HT9GUA 3, and by the European Commission through the Research 
Training Network Quantum Chromodynamics and the Deep Structure of Elementary Par- 
ticles under Contract No. ERBFMRXCT980194. The work of A.A.P. was supported in 
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Note added 

After including the 0(a^/3) contribution from Eq. (|l^), which was missed in the pre- 
vious version of this letter, we find agreement with the analytical results for the p-Ps and 



o-Ps decay widths of Ref. |26| and with the numerical result for the o-Ps decay width of 



Ref. 1^ . However, this correction is of insignificant size and immaterial for our numerical 



estimates of the Ps decay widths. 
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